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ABSTRACT 



O 

We here propose a two-step model for gamma-ray bursts (GRBs) associated with su- 
pernovae. In the first step, the core collapse of a star with mass > 19M© leads to a 



\ massive neutron star and a normal supernova, and subsequently hypercritical accretion 
of the neutron star from the supernova ejecta may give rise to a jet through neutrino 

^ \ annihilation and/or Poynting flux along the stellar rotation axis. However, because of 

OO 
O 

o> 
o> 



too much surrounding matter, this jet rapidly enters a nonrelativistic phase and evolves 
to a large bubble. In the second step, the neutron star promptly implodes to a rapidly 



^ ■ rotating black hole surrounded by a torus when the mass of the star increases to the 
maximum mass and meanwhile its rotation frequency increases to the upper limit due 



to the accreted angular momentum. The gravitational binding energy of the torus may 
be dissipated by a magnetized relativistic wind, which may then be absorbed by the 
supernova ejecta, thus producing an energetic hypernova. The rotational energy of the 
black hole may be extracted by the Blandford-Znajek's mechanism, leading to another jet. 
This jet is relatively free of baryons and thus may be accelerated to an ultrarelativistic 
phase because the first jet has pushed out of its front matter and left a baryon-free exit. 
We expect that the second jet generates a GRB and its afterglow. Our two-step model 
may alleviate the baryon-contamination problem suffered possibly from in the hypernova 
models. Furthermore, this model not only accounts for association of several GRBs with 
supernovae but also explains well the features of the afterglows of these bursts. 

Subject headings: gamma-ray: bursts - supernovae: general - stars: neutron - black holes: 
physics 
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1. INTRODUCTION 



It has been widely believed that gamma-ray burst (GRB) events do indeed occur at 
cosmological distances (Metzger et al. 1997; Kulkarni et al. 1998a, 1999; Andersen et 
al. 1999), which implies that a successful model for progenitors of cosmological GRBs 
must satisfy two essential requirements: (1) The model should produce an extremely 
relativistic fireball, which should subsequently emit an amount of gamma-ray isotropic 
energy E y ~ 10 51 -10 54 ergs implied by the observed fluences and the cosmological distance 
scale. The recent multi-wavelength observations of GRB afterglows support the so-called 
fireball shock model (Piran 1999). According to this scenario, GRBs are produced as a 
result of internal shocks when fast moving shells catch up with slower shells that were 
ejected at earlier times. The relative kinetic energy of motion of the shells is converted into 
observed gamma-ray emission in relativistic shocks via synchrotron radiation or inverse- 
Compton scattering mechanisms. The observed afterglow emission is produced when the 
shell decelerates as a result of interaction with the ambient matter. Based on the fireball 
shock model, several authors (Panaitescu, Spada & Meszaros 1999; Kumar 1999; Lazzati, 
Ghisellini & Celotti 1999) recently found that the efficiency for producing gamma-rays 
in internal shocks is a few percent. Therefore, the isotropic energy of fireballs in some 
bursts (e.g., GRB 990123) must be up to a few times 10 55 ergs (Kulkarni et al. 1999). 
If anisotropic emission with a beaming factor of Ail/An ~ 0.01 is assumed, this energy 
can be reduced to Ej Ct ~ a few times 10 53 ergs. (2) The rapid variability of GRBs and 
their nonthermal spectra (Woods & Loeb 1995) and the low radiative efficiency in internal 
shocks (Lazzati et al. 1999) requires that the Lorentz factor of the fireball be 100-1000. 
This implies that the fraction of contaminating baryons must be less than 1%. If the 
emission is anisotropic, the mass of loading matter AM < 0.01£'j et /c 2 ~ 1O~ 3 M . 

Two currently popular models of GRB progenitors are the merger of two compact ob- 
jects (neutron stars and black holes) and the collapse of massive stars. The former model 
is the plausibly baryon-clean one, but compact objects would be expected to have such 
a significant space velocity that their merger would take place outside their birthplaces 
(Paczyhski 1998; Bloom, Sigurdsson & Pols 1999). The observational evidence for the as- 
sociation of several GRBs with star forming regions then provided weak evidence against 
the compact object merger as GRB progenitors and favored massive star progenitors. The 
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population synthesis of Fryer, Woosley & Hartmann (1999) supported this conclusion. 

The massive star progenitor model has become more favorable since the discovery of 
the Type Ic supernova SN 1998bw in the error box of GRB 980425 (Galama et al. 1998). 
The high energy inferred for the optical supernova, (2 — 3) x 10 52 ergs (Iwamoto et al. 
1998; Woosley, Eastman & Schmidt 1999), and the high expansion velocity inferred for the 
radio supernova (Kulkarni et al. 1998b) strengthen the GRB-SN connection. Recently, 
this connection was further confirmed by Bloom et al. (1999), Reichart (1999) and Galama 
et al. (1999), who found the dramatical brightening and extreme reddening of the optical 
afterglows of GRB 980326 and GRB 970228 at late times, respectively. However, it is quite 
difficult to understand such a connection. From the above second essential requirement, 
too many baryons possibly exist in the vicinity of the collapsing core in the hypernova 
models (Woosley 1993; Paczyhski 1998) so that an ultrarelativistic jet forming during 
the collapse of the core rapidly become nonrelativistic. This conclusion is consistent with 
the numerical studies (MacFadyen & Woosley 1999): an ultrarelativistic jet converts to 
a nonrelativistic large bubble. A simple reason for this may be that a large amount of 
radiative energy is released impulsively but the mass of contaminating baryons is in fact 
of the order of M,(Afi/4vr) ~ O.lM (M*/lOM )(lO 2 Afi/47r), where M* is the mass of 
the matter around the core. 

In this Letter, we propose a scenario, in which a supernova explosion first may produce 
a massive neutron star, and about two hours later the star will start to accrete the fall-back 
supernova ejecta at a hypercritical rate, while an accretion disk will form near the stellar 
surface due to its large angular momentum. In particular neutrino annihilation and/or 
Poynting flux along the rotation axis of the star may lead to a jet which will pushes 
out of its front matter. The mass of the accreting neutron star will eventually reach the 
maximum mass about several hours after the supernova explosion and thus will promptly 
collapse to a rapidly rotating black hole surrounded by a torus. The gravitational binding 
energy of the torus, which is of the order of several 10 52 ergs, may be dissipated into the 
supernova ejecta, which may in turn give rise to an energetic hypernova. Another jet with 
energy of a few times 10 53 ergs and with low-mass baryon contamination will occur along 
the rotation axis of the black hole by extracting its rotational energy via the Blandford- 
Znajek's (1977) mechanism. We expext that the second jet can produce a GRB and its 
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afterglow. 

Delayed formation of black holes in supernovae has been widely discussed (Woosley 
1988; Chevalier 1989; Brown & Weingartner 1994; Brown & Bethe 1994; Woosley & 
Weaver 1995; Fryer 1999). Recently, it has been shown numerically that the accretion 
of such a black hole may lead to a relativistic jet required by a typical GRB (Woosley, 
MacFadyen & Heger 1999; MacFadyen, Woosley & Heger 1999). The present model is 
an analytical one, in which we suggest that hypercritical accretion of a newborn neutron 
star and delayed formation of a black hole could produce two discrete jets. The first jet 
will push out of its front matter and leave an exit for the second jet, which could thus be 
relatively free of baryons. 

2. HYPERCRITICAL ACCRETION OF A NEWBORN NEUTRON STAR AND 

FORMATION OF A BUBBLE 

It is well known that the core collapse of massive stars with 1O-25M produces neutron 
stars accompanying Type II supernovae. Timmes et al. (1996) numerically studied the 
initial mass function of newborn neutron stars and found that their initial mass distribu- 
tion is bimodal with peaks at 1.27 and 1.76M . The principal reason for this bimodal 
distribution is the differnce in the presupernova structure of stars above and below 19M , 
the mass seperating stars that burn carbon convectively from those that produce less 
carbon and burn radiatively. Here we consider neutron stars with 1.76M Q as the startimg 
point of our work. When such a massive neutron star first forms in a supernova explosion, 
it is surrounded by a dense gas (supernova ejecta), some of which falls onto the neutron 
star and cools by neutrino emission (Colgate 1971; Zeldovich, Ivanova & Nadezhin 1971). 
It is the neutrino emission that allows accretion of the star at a high rate. From sim- 
ple analytical arguments, Chevalier (1989) and Brown & Weingartner (1994) estimated 
a lower limit to steady neutron star accretion with neutrino losses assuming spherical 
symmetry. However, the accreted matter may have a large angular momentum which 
leads to an accretion disk. In this case, the lower limit with neutrino losses is estimated 
as M cr ~ 1.1 x lO^Moyr- 1 (Chevalier 1996). 

The supernova explosion scenarios involve an outgoing shock wave. When this shock 
enters the hydrogen envelope, the deceleration of matter occurs. This deceleration sharp- 
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ens into a reverse shock. When the reverse shock reaches the neutron star surface, 
the star starts to accrete the fall-back supernova ejecta at a hypercritical rate. The 
time scale for the reverse shock to pass through the core is about the same as the 
time that the outgoing shock front takes to reach the stellar surface, which is t ~ 
2hi(R*/3 x lO 12 cm)(M ej /lOM ) 1 / 2 (£; sn /lO 51 ergs)- 1 / 2 , where R* is the presupernova stel- 
lar radius, M e j is the mass of the ejected matter and E sn is the supernova explosive energy 
(Shigeyama, Nomoto & Hashimoto 1988). For SN1987A, t is about 2 hr (Shigeyama et 
al. 1988). In order to derive the rate of mass accretion onto the neutron star after to, 
we follow Brown & Weingartner (1994). By assuming the neutron star to be at rest with 
respect to its ambient matter and using the Bondi's (1955) spherical accretion theory, we 
obtain the accretion rate 

M = 5 . 63 x 10- (")' (j^Y^ r'Meyr-, (1) 
\1.8Mq) V10 8 cms-V \1jt J OJ ' w 

where M is the neutron star mass (~ 1.8M Q ) and Vf is the final velocity (after being slowed 
down by the reverse shock) of the carbon-oxygen core. For SN1987A, v / ~ 6 x 10 7 cms _1 
(Woosley 1988). In deriving equation (1), we have assumed that the ambient matter of the 
neutron star is radiation-dominated due to the effect of the outgoing shock. Owing to this 
assumption, our accretion rate is slightly larger than that of Brown & Weingartner (1994). 
The time at which radiation significantly affects accretion can be estimated based on 
M(t CT ) = M cr . Hence, we obtain this timescale t CT = 0.7vjlyr, where Vf$ = t>//10 8 cms _1 . 
Below, we take i>/,8 = 1- When t < t cr , the accretion is hypercritical and the total accreted 
baryon mass is given by 

AM acc = /* Mdt = 1.81M ^ 5/8 [(t o /lhr)- 7 / 8 - (t/lhr)- 7 / 8 ]. (2) 

J t •" 

For the modern realistic equation of state for neutron matter chosen in the next section, 
an accreting neutron star with initial mass of 1.76M will collapse to a rapidly rotating 
black hole when AM acc = O.55M . 

Before the collapse, the accreted matter forms a disk near the neutron star because the 
accreted angular momentum may be up to 3 1 ^ 2 R s c per gram, where R s is the Schwarzchild 
radius (Woosley & Chevalier 1989). The temperature of the accretion disk near the 
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neutron star can be estimated by the following equation: r/Mc 2 = 



e u An d Rl s )/A, where 



r\ is the efficiency for the conversion of the gravitational energy to heat (~ 0.1), R ns is 
the neutron star radius (~ 10 6 cm), Af2 d is the solid angle of the disk, and e v = 1.0 x 
10 25 (T/MeV) 9 ergs _1 cm" 3 is the neutrino pair energy production rate per unit volume 
(Dicus 1972). For typical values of these parameters (e.g., A£l d ~ 3), we obtain T ~ 6 
MeV. The total energy for neutrino losses is approximated by E v = GMAM acc /i? ns ~ 
3 x 10 53 ergs. Since anisotropic neutrino emission takes place due to the effect of the disk 
(Kluzniak 1998), neutrino annihilation along the rotation axis of the neutron star leads 
to a jet. By using the efficiency of neutrino annihilation, x ~ 0.3% (Goodman, Dar & 
Nussinov 1987; Kluzniak 1998), we obtain the energy of the jet 



It should be pointed out that this efficiency is much larger than that of Popham, Woosley 
& Fryer (1999), in which an accretion disk surrounding a black hole is advection-dominated. 
In the present case, however, matter is accreted onto the neutron star surface where the 
conversion of the gravitational energy to heat (77) is much more efficient. Since this ther- 
mal energy is released via neutrino emission, the efficiency of neutrino annihilation should 
increase substantially. In addition to the neutrino annihilation mechanism, there is an- 
other possible mechanism to produce a jet proposed by Katz (1997), in which the magnetic 
field amplified by the differential rotation of the disk may result in a strong Poynting flux. 
The jet produced by these mechanisms will push its front baryonic matter whose veloc- 
ity is given by u jet = (2£ jctil /AM') 1/2 ~3x 10 9 (M*/10M o )- 1 / 2 (l0 2 Afi/4 7 r)- 1 / 2 cms" 1 , 
where AM' = M*(Afi / 'Att) is the mass of baryons loading with the jet. It should be 
noted that this velocity is much larger than that of the outgoing shock. As numeri- 
cally studied by MacFadyen & Woosley (1999), this nonrelativistic jet will expand lat- 
erally, producing a bubble. When t ~ 5 hr, at which the collapse of the neutron star 
may take place as argued in the next section, the radius of the bubble is estimated as 
Rb ~ Ujet(*-*o) ~ 3 x 10 13 (M,/10M o )" 1 / 2 (10 2 Afi/47r)- 1 / 2 cm, which is much larger than 
the presupernova stellar radius R*. 

3. COLLAPSE OF THE NEUTRON STAR AND GENERATION OF A GRB 



'jet,l 



E v \ ~ 10 5i ergs. 



(3) 
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It has been argued by Brown & Bethe (1994) and Bethe & Brown (1995) that hypercrit- 
ical accretion of neutron stars in supernovae (e.g., SN1987A) may lead to collapse of the 
stars to low-mass black holes. These authors considered a soft equation of state like kaon 
condensation for neutron matter as the starting point of their discussions. Here we adopt 
a more realistic modern equation of state named UV14+TNI (Wiringa, Fiks & Fabrocini 

1988) . This equation of state gives the property of a nonrotating neutron star at the max- 
imum mass: maximum gravitational mass M max = 1.84M , corresponding baryon mass 
Mf, = 2.17M Q , radius _R ns = 9.51 km and moment of initia / = 1.5 x 10 45 g cm 2 (Wiringa et 
al. 1988). According to this property, the maximum frequency of rotation for this neutron 
star equation of state is given by Vt max = 7.84 x lO 3 (M max /M ) 1 / 2 (i^/lOkm) -3 / 2 s -1 = 
1.17 x 10 4 s _1 (Cook, Shapiro & Teukolsky 1994). The property of a rotating neutron 
star at this maximum frequency and at the allowable maximum mass is as follows (Cook 
et al. 1994): maximum gravitational mass M max = 2.19M , corresponding baryon mass 
M b = 2.55M , radius R ns = 12.7 km and angular momentum J = 2.85 x 10 49 ergs or the 
Kerr rotation parameter a = cJ/GM^ ax = 0.67. 

After having these properties, we discuss implications of hypercritical accretion of a 
newborn neutron star with gravitational mass of 1.76M . First, hypercritical accretion 
may produce and maintain a large bubble, as argued in the last section. 

Second, hypercritically accreted matter may rapidly submerge the magnetic field of 
the neutron star (Muslimov & Page 1995). Once the accreted mass reaches O.O1M , the 
buried magnetic field ohmically diffuses out after ~ 10 s yr (Geppert, Page & Zannias 
1999). This implies that the magnetic field of the neutron star could be always weak in 
the accretion timescale. 

Third, hypercritical accretion may spin up the neutron star (Woosley & Chevalier 

1989) . A 1.76M neutron star has the baryon mass of 2.OM and thus this star needs 
to accrete matter with mass of O.55M to become a maximum rotating neutron star 
with maximum mass. As suggested by Woosley & Chevalier (1989), the accreted angular 
momentum from the mixed mantle and helium core of the ejecta may be as large as 3 1//2 i? s c 
per gram corresponding to 5.54 x 10 48 ergs per O.1M accreted mass. Thus, after accreting 
O.55M mass, the neutron star can obtain the angular momentum of 3 x 10 49 ergs. But, 
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the maximum angular momentum of a rapidly rotating neutron star at maximum mass 
is only 2.85 x 10 49 ergs. The remaining angular momentum is A J = 0.15 x 10 49 ergs. 
How is AJ dissipated? Fortunately this remaining accreted angular momentum can be 
easily carried away by gravitational radiation (Wagoner 1984). Owing to hypercritical 
accretion, therefore, the mass of the neutron star not only reaches the maximum mass 
but it is spun up to the maximum rotation frequency. 

Finally, once the mass of the accreted matter reaches 0.55M Q , the neutron star will 
promptly collapse to a rapidly rotating black hole. From equation (2), we can estimate 
the time at which the collapse will occur: t ~ 5.1 hr. The resulting black hole has the 
rotational energy 



where f{a) = 1 - J (I + y/1 - a 2 )/2 pa 0.067 at a pa 0.67, and the mass of the black hole 



has been taken to be 2M . After the collapse, not all mass will be immediately accreted; 
the outmost layers with a small fraction (a few percent) of the total mass, in fact, have 
centrifugal accelerations beyong the local gravitational attraction (Vietri & Stella 1998), 
leading to a torus. In addition, a small amount of accreted ejecta still stagnate in the torus. 
Thus, the total mass of the torus may be M t > O.1M . The presence of the torus will give 
rise to the following two effects: (i) The binding energy of the torus is Ef, = GM-QuMt/ Rt = 
3x 10 52 eYgs(M t /0.1M & )(R t /20 km) -1 , where Rt is the typical radius of the torus. Because 
the magnetic field of the torus may be amplified to 10 15 G due to differential rotation in 
the torus, this binding energy can be dissipated into a magnetized relativistic wind (Usov 
1994; Meszaros & Rees 1997; Katz 1997). Because such a wind may be easily absorbed 
by the outgoing supernova ejecta as argued by Dai & Lu (1998a, 1998b), almost all of 
the binding energy may convert to the expansion energy of the ejecta, which may further 
give rise to a supernova with a much brighter optical luminosity, (ii) In the presence of 
the torus, the rotational energy of the hole can be extracted by the Blandford-Znajek 
(1977) mechanism (Meszaros, Rees & Wijers 1998). The power for this mechanism is 



Pbz = 1.7xlO 51 a 2 /(a)(M BH /M ) 2 ( J B/3xlO 15 G) 2 ergs- 1 = 2.5x 10 50 (5/3x 10 15 G) 2 erg S - 1 



(Lee, Wijers & Brown 1999). The rotational energy of the black hole will be dissipated 
in the timescale of ~ 10 3 s if the magnetic field strength is of the order of 3 x 10 15 G. 



E rot = f(a)M BH c 2 pa 3 x 10 53 ergs, 



(4) 
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We believe that this mechanism will produce an ultrarelativistic jet because the mass of 
matter loading with this jet can be estimated by 



AM = AM'(Aft/47r) = lO~ 3 M (M,/lOM )(lO 2 Afi/47r) 2 . 



(5) 



In fact, since the angular momentum of the presupernova star may be rather large, the 
centrifugal force should reduce the density of matter along the rotation axis and thus the 
above estimate is an upper limit. After a short acceleration phase, the Lorentz factor of 
the jet should be 



which is consistent with the constraints from the rapid variability of the light curves and 
the nonthermal spectra of GRBs (Woods & Loeb 1995). The collisions among the shells in 
such a jet will produce a GRB and subsequently the deceleration of the jet in its ambient 
medium will result in an afterglow. 



We have proposed a two-step model for GRBs associated with supernovae. A hyper- 
nova, a much brighter supernova, may occur in our model. However, this model is clearly 
different from the current hypernova models (Woosley 1993; Paczyhski 1998; MacFadyen 
& Woosley 1999). In the latter models, the core of a massive star directly collapses to a 
black hole. Recently, Cen (1998) and Wang & Wheeler (1998) proposed models for GRB- 
SN association, in which the matter above the neutrinosphere in a small cone around some 
special axis of a newborn neutron star is assumed to be preferentially first blown out of 
the deep gravitational potential well of the star in order to avoid too many baryons con- 
taminating a subsequently resulting jet. Our model may provide a plausible way of how 
such an empty cone is produced: neutrinos from the hypercritical accretion disk annihi- 
late to electron/positron pairs which form the first jet to push its front baryons and leave 
an exit for the second jet. Therefore, our model may avelliate the baryon contamination 
problem. 



r = E Iot /{AMc 2 ) > 15O(M,/lOM )- 1 (lO 2 Afi/47r) 



-2 



(6) 



4. DISCUSSIONS 
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In our model, GRBs are naturally associated with supernovae because the former 
phenomenon takes place several hours after the latter phenomenon. This is consistent 
with the analysis of Iwamoto et al. (1998), in which the time of core collapse coincides 
with that of SN1998bw within (+0.7, —2.0) days. The unusually large explosive energy 
of SN1998bw, (2 -3) x 10 52 ergs (Iwamoto et al. 1998; Woosley et al. 1999), is very close 
to the gravitational binding energy between the black hole and torus, ~ 3 x 10 52 ergs. 
This consistency seems to support our model. 

The features of some afterglows can also be explained in the context of our model. For 
example, the optical afterglow of GRB 980326 rapidly decayed as oc t- 2 0±01 in the first 
two days, subsequently brightened dramatically and reddened significantly, and finally 
declined (Bloom et al. 1999). These features were also seen in the light curve of the 
afterglow of GRB 970228 (Reichart 1999; Galama et al. 1999). The late-time afterglows 
were widely believed to be the contribution of unusually brighter supernovae, while the 
rapid decay of the early afterglows was analytically argued to be due to sideways expansion 
of jets (Rhoads 1999; Sari, Piran & Halpern 1999). Another interpretation for the rapidly 
decaying afterglows may be that a relativistic shock expanding in a dense medium has 
evolved to a nonrelativistic phase (Dai & Lu 1999a, 1999b). In our two-step model, this 
dense medium may be a presupernova steller wind. 

In summary, our two-step model may alleviate the baryon-contamination problem suf- 
fered possibly from in the hypernova models. This model not only accounts for association 
of several GRBs with supernovae but also explains well the features of the afterglows of 
these bursts. 

We would like to thank J. S. Bloom, B. Hansen, A. MacFadyen and S. E. Woosley for 
their comments and discussions. This work was supported by a RGC grant of Hong Kong 
government and the National Natural Science Foundation of China. 
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